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Abstract. We extend earlier results on conditioning of super-Brownian motion to gen- 
eral branching rules. We obtain representations of the conditioned process, both as an 
/j-transform, and as an unconditioned superprocess with immigration along a branching 
tree. Unlike the finite-variance branching setting, these trees are no longer binary, and 
strictly positive mass can be created at branch points. This construction is singular in 
the case of stable branching. We analyze this singularity first by approaching the stable 
branching function via analytic approximations. In this context the singularity of the 
stable case can be attributed to blowup of the mass created at the first branch of the 
tree. Other ways of approaching the stable case yield a branching tree that is different 
in law. To explain this anomaly we construct a family of martingales whose backbones 
have multiple limit laws. 



1. Introduction 

The ip- super-Brownian motion X t is a measure valued diffusion with branching mecha- 
nism given by 

POO 

ip(X) = a x X + a 2 A 2 + / [e" Ar — I + Ar] ir(dr), 



where A > 0, a x e R, a 2 > 0, and 7r(-) is the associated Levy measure. More precisely, 
the log-laplace functional of X t , 



u(t,x) = -log Eg. exp{~ j (f>(y)dXt(y)) 



solves the initial value problem 

°- = -^u-^[u\ u(0, •) = #•), 

whenever is a non-negative bounded continuous function. These diffusions can be 
realised as scaling limits of a system of particles that perform branching Brownian motions. 
The exit measure X D of such a process from a bounded domain D is a random measure 
on the boundary of D, supported on the set of points where the particles first exit the 
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domain. It is known that for d > 3 that the random measure does not see points (see Le 
Gall [2Q] and Dynkin [5]). 

In Salisbury and Verzani [29] and [30] the exit measure of super Brownian motion, X D with 
critical binary branching (that is, with branching mechanism = 2A 2 ) is conditioned to 
charge small balls Af — B(zi,e) D dD, where Z\, . . . , z n G dD. Letting e — >■ they obtain 
a conditioned process, which is a martingale transform of super Brownian motion by a 
"polynomial" martingale of degree n. They describe this process in terms of a "backbone" 
consisting of a binary tree that realizes the trajectory of the mass that reaches z\, . . . , z n . 
These results do not generalize to a stable branching mechanism "0(A) = cA 1+/3 , since 
"polynomials" of the exit measure will not even be integrable, let alone give martingales, 
when n > 2. Understanding why formed the primary motivation for this paper. 

Our aim is to understand conditioning based on general branching mechanisms well 
enough to analyze how these martingales blow up as we approach the stable case. Consider 
a bounded domain D C M d , when d > 3. Let e > and fix Zj G dD for i — 1, 2, . . . , n. As 
above, define Af to be a ball on dD of radius e. We condition ^-super Brownian motion 
to hit balls Af and obtain both the martingale transform that represents this process, 
and the probabilistic representation of this process in terms of immigrating mass along 
a branching "backbone" . We explicitly describe the evolution of the backbone tree, the 
manner mass is generated along the backbone, and the way it evolves afterwards (See 
Theorem 13.11 in Section [3]). Unlike the results of [29] and [30], we now have to handle 
the probabilities of multiple branches and the distribution of positive mass created at the 
branch points of the tree. 

We first take the limit as e — > 0, when d > 4 and ip is a real-analytic function. We 
establish that the limit exists and is a martingale change of measure (see Theorem 14.21 
and Theorem 14.11 in Section 14.21 ). The proof here follows the road map laid out in [2D] 
but the significant estimates appear to require more delicate arguments (Lemma 14.21 and 
Lemma I4~3l) . The limiting backbone produced is a tree with precisely n leaves (see Section 
14. 3p . but is no longer binary. Next we let the branching mechanism approach the stable 
case (V'(A) = cA 1+ ^). With this particular order of taking the two limits, the backbone 
remains well behaved, as does the mass creation at non-branch points of the backbone. 
However, the mass created at branch points gets arbitrarily large, resulting in an explosion 
that allows us to pinpoint the source of the blowup (See Section I4~4l and Theorem |4.4p . 
We also consider other ways of approaching the stable case, for example, to simply let 
e — > with branching mechanism "0(A) = cA 1+/3 . In this case, the backbone remains well- 
behaved but has a different law than in the previous mechanism. To explain this anomaly 
we construct a simpler family of martingales whose backbones have multiple limit laws, 
interpolating between analogues of both types obtained above (See Section |4~5I) . 

The decomposition of the conditioned superprocesses in terms of an "immortal backbone" 
has been considered in the literature, in other contexts. In particular, /i-transforms of 
critical super-Brownian motion have been studied by Roelly-Coppoletta and Roualt [27j . 
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Evans and Perkins [17], and Overbeck [21], [22], while the immortal particle representation 
was originally discovered in Evans JT5J. Other studies include Serlet [21] and Etheridge 
[H]. Salisbury and Sezer [28] and Verzani [33] consider more general conditionings for 
binary branching. Moras [23] considers specific classes of unbounded domains D, again 
for binary branching. 

Versions for non-binary branching were studied in Etheridge and Williams [15] and Kypri- 
anou et al [19] (we learned of the latter after completing our research). In [15], the super 
Brownian motion on all of M d , with stable (1 + /3) branching, is conditioned on survival 
until some fixed time T. This backbone has a Poisson number of immortal trees (con- 
ditioned on there being at least one), along which mass (conditioned to die before time 
T ) is immigrated. The rate of immigration is random and there is additional immigra- 
tion whenever the immortal tree branches. In the limit as T — > oo , the immortal trees 
degenerate to the Evans immortal particle and the immigration (of unconditioned mass) 
along the particle is dictated by a stable subordinator. In [19], to study the travelling 
wave equation associated to the parabolic semi-group equation of ^-super-Brownian mo- 
tion, the authors show a similar immortal backbone on all of R (which they call a "spine 
decomposition"). 

1.1. Layout of the paper. The rest of the paper is organised as follows. In the next 
section we discuss preliminaries with regard to conditioned diffusions, ^-super Brownian 
motion and potentials. The Palm formula in Lemma 12.31 and identity for potentials 
in Lemma 12. 5\ presented here, are used significantly in the proof of main results. In 
Section [3] we prove Theorem 13.11 In particular, we define a martingale change of measure 
representing various conditionings and for each provide a description of the associated 
branching backbone representation. 

In Section HI we consider a specific condition namely that of the exit measure charging 
finitely many points on the boundary. The work here (as explained earlier) is divided 
into three parts. First for ^-analytic, we condition the exit measure to hit balls of radius 
e and establish a limit as e — > 0. This is proved in Theorem 14.11 and Theorem 14.21 A 
limiting backbone is then described in Section 14.31 Secondly, in Section 14.41 we consider 
the second limiting procedure of ^-analytic to approximate i(j(X) = c\ 1+/3 for < (3 < 1. 
We explain the explosion effect precisely. Finally in Section I4.5[ we begin by explaining 
other possible limits (if the order of limits done earlier are interchanged). We conclude 
by constructing a family of related martingales whose backbones have multiple limit laws. 
These interpolate between analogues of the two types of limits obtained earlier. 

1.2. Acknowledgements. Research done in this paper was supported in part by NSERC. 
During the completion of this work: the authors visited the University of British Columbia; 
the Fields Institute; Siva Athreya visited York University in Toronto; and Tom Salisbury 
visited the Indian Statistical Institute, Bangalore. We would like to thank all these places 
for their kind hospitality. 
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2. Preliminaries 



2.1. Notation. For each set A, let \A\ denote its cardinality, and let V(A) denote the 
collection of partitions of A. Impose an order on A. Then for any a G V(A), we may 
order the sets in a by their smallest elements. Letting a(j) denote the jth element of o 
in this order, we may switch at will between the following two notations: 



U(X D ,v c )=H(X D ,v^). 

C&a j=l 



The following elementary combinatorial result will be convenient 
Lemma 2.1. Let A C B C C. be subsets of {1,2 ... ,n} . Then 

-1)1*1 = (-1)'°'U, C 



E ( 

AcBcC 



Proof. Both sides equal (—1)^1(1 — l)'*"^'. See Lemma 2.1 of 



□ 



2.2. Facts about conditioned diffusions. First we recall some familiar formulae for 
conditioned Brownian motion. 

Let B be rf-dimensional Brownian motion started from x, under a probability measure 
P x . Write Tr> = Td(5) for the first exit time of B from D. Let g : D — > [0, oo) be bounded 
on compact subsets of D, and set 




Let £t be a process which, under a probability law P£, has the law of a diffusion with 
generator L g started at x and killed upon leaving D. In other words, £ is a Brownian 
motion in D, killed at rate g. Write ( for the lifetime of £. Then 

(2.1) P*(£ t eA 7 C>t)=P x (exp- J g(B,)ds,B t eA,T D >ty 

Let U 9 f(x) = J °° P£(f(£,t)1-{c>t}) dt be the potential operator for L g . If g = we write U 
for U 9 . If < u is L 9 -superharmonic, then the law of the w-transform of £ is determined 
by the formula 

for $(£) G cx{£<,; s < t}. Assuming that < u < oo on D, this defines a diffusion on D. 
If it is L 9 -harmonic, then it dies only upon reaching 3D. In fact, the generator of the 
^-transform is 

L g . u f = -L g (uf) = l-Af + -\7u ■ V/. 

u 2 u 
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If u = U 9 f for some / > (that is, if u is a potential) then the u-transform dies in the 
interior of D, and P^' u satisfies 

prmo) = 4t r P2m<t)m)Mc>t})dt, 

u \ x ) Jo 

where £< 4 is the process £ killed at time t. See [2]. 

More generally, if L g u = —f then by breaking u into a potential plus a harmonic function 
one has that 

(2.2) prm)h«r D} ) = 4t f°° p ^<t)f&)M<>t}) dt. 

u \ x ) Jo 

2.3. i/i super Brownian-motion. Let 

P oo 

(2.3) = a 1 X + a 2 X 2 + / [e~ Ar - 1 + Ar] ?r(dr), 

where A > 0, a 2 > 0, and 7r(-) is the associated Levy measure. We will assume the 
following: 

POO 

(2.4) ai > 0; / min(r, r 2 ) n(dr) < oo. 

Levy process exist under the weaker condition J °° min(l, r 2 ) n(dr) < oo, and indeed 
there is a well known construction of continuous state branching processes (CSBP) as 
time-changes of Levy processes. The stronger moment condition assumed above can be 
thought of as a condition for this CSBP (or equivalently, this time change) not to blow 
up in finite time. See [18]. It is easily seen that ip^(X) = 4^^(X) exists for all n and 
A > 0. 

Let D be a domain in M. d . Take £ t to denote Brownian motion on D, and let X t be the 
^-super Brownian motion on M. d . Let X D be the associated exit measure on 3D. In [12], 
Dynkin constructs this object, assuming a± = 0. But by Dawson's Girsanov theorem (see 
section 10.1.2 of [7]), applying Dynkin's construction to Brownian motion killed at rate 
ai yields precisely the desired ^-super Brownian exit measure. Note that we are using 
the condition ai > at this point. For a\ < and tx = one could in fact produce 
a superprocess that survives forever with positive probability, by taking D large enough 
(see [13]). Thus the exit measure could be infinite in that case. 

N x will denote the excursion measure. LeGall used this measure extensively in the case 
"0(A) = 2A 2 (see [21]). For the general case see Dynkin [12]. A key benefit to working under 
Rj, is that genealogies simplify - all mass descends from a single massless initial individual. 
The price one pays for this simplification is that N x is an infinite measure. Only events 
that involve extinction at short times receive infinite mass, so N X (X D ^ 0) < oo. It can be 
realized by having the superprocess start with initial value jd x under a probability measure 
P 7< 5 X , sending 7 J, 0, and renormalizing F^s x to obtain a non-trivial limit. Alternatively, 
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the superprocess with initial value \x (under a probability P M ) can be realized in terms of 
a Poisson random measure having J N x fi(dx) as intensity. See section [3] 

Let = exp —(X D , 0). When convenient we will also denote this e D (0) or just e(0). Let 
Dk be a sequence of smooth subdomains increasing to D and denote X Dk by X k and e® k 
by e$. Set 

(2.5) M t {f) = exp ( - jf <0'(%(1 - /)) ds] 
Since ai > we have > when A > 0, so in particular A/*(/) < 1 for / > 0. 
We will need the following results. 

Lemma 2.2. Assume f!2.4[) and /ei D be a domain in ~R d . 

(a) Lei T C <9L>. TYien o(y) = f\(X D (r) > 0) sato/ies \Ag = ip(g). 

(b) Lei g be a non-negative solution to \Ag = tp(g), and let Dk be an increasing 
sequence of smooth subdomains of D. Then for each k, 

N x (l-e k g )=g(x). 

(c) X x ({X D ,<j>)ef) = E x m TD )K D (ef)) for f, > 0. 

Proof. Follows from Theorem 4.2.1 in [10], Theorem 1.1. in Chapter 4 of [12] and Theorem 
11.7.1 in 0. 

□ 

For m > 2 an integer, ?/ G R d , and measurable we define 

(2.6) b(m, 0, y) = (-l)"V (m) (N w (l - e,)) . 
Note that for A > 

{ai + 2a 2 A + J °° r[l - e" Ar ] ir(dr), m = 1 
2a 2 + / OO r 2 e- Ar 7r(rfr), m = 2 

(-l) m J °° r m e" Ar vr(dr), m > 3, 

so in particular, 6(m, 0, w) is well defined and non-negative, for m > 2 and for any > 
such that N y ((X D , 0} > 0) > 0. It is decreasing in 0. 

We will need the following assumption on ip for certain results, which, among other things, 
makes the latter qualification unnecessary. 

(AI) 3A > such that f™ e rX ° vr(dr) < oo. 
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Under (Al) we have that J °° r n n(dr) < oo for n > 2, and by dominated convergence, 



'-* J \j POO 

(2.8) iP(\) = ai \ + a 2 \ 2 + Y(-iy— r^{dr) 



for A < An. Thus in this case, (12.71) will hold as well for A = 0. 



Recursive Palm formulae for moments have a long history, and in this context are due 
to Dynkin - see Theorem 1.1 of [12] or Theorem 4.1 of [11]. We work instead with a 
formulation along the lines of Lemma 2.6 in [29] . 

Lemma 2.3. Assume (12. 4p and let N = {1, 2, 3, . . . n}, n > 2. Let D C W 1 be a domain, 
(j) > 0, and let £ be a Brownian motion in D with exit time r. Let {vi} be a family of 
positive measurable functions. Then 



N x (e ( ,l[(X D ,v l ))=E x ( I ^t(e 4> )b(\f3\, ( f ) ^ t )l[N^l[(X D ,v l ))dt 

\B\>2 



Proof. Let N* = {2, 3, . . . , N}. 

N x (e^]J(X D 1 v l })=N x ((X D 1 v 1 }e^l[(X D ,v l }) 

_ x d d d 



:-iy 



d\ 2 d\ 3 d\ n 



N x ((X D ,v 1 )e( ( j ) + J2^i))^ 

A 2 = ...=A n =0 j = 2 



where we differentiate under the integral sign using monotone convergence. Using (c) of 
Lemma I2T21 i.e. the one-dimensional case of the Palm formula, the above expression equals 



n 

£zM£r)X(e(0 + ^ A ^))) 
A 2 =-=A n =0 j = 2 
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Differentiating this expression under the integral sign (given the conditions on n this can 
be easily justified) and using the definition of b, this equals 



(2.9) 



S,UW(e*) E 11/ E (-l) lftl+ V lftl+1 N6,(l 



- x 



Ae/3j i€A 



Standard integration manipulations then show that the above equals 



e x ( v 1 (t T )K(e i> ) e e f e n N ft fc (^n< xD ' t, «» : 



n ( f e + x > & ^) n (e* n< xD »<» rft 



e A e e / r e n N ^^iK xD '^) : 

\<re7>(JV*) fc=l ^ k eV(<r(k)) A€/3 k i<=A 

XE X (vi(er)^r(e^n f E KI&I + I^&JX 



Ag/3,- i£A 
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where denote the filtration of £ tk . Applying the Markov Property at time tf., this equals 

\a£V(N*) fc=l /3 k eV(a(k)) A€/3 k i£A 

x n^^n^'^)^) 

Setting M = <r(/c) and summing over N* this becomes 



^( e [mm e km + 1,0,6) n%(^n^'^)) x 



x 

pep{N*\M) \ j=i" u ftePMj)) 

x II %(e«n^' u «)) ds ) dt 
Using the identity obtained in (12.91) . this equals 

eJ r mm y Ki/3i+i,0,6)n^^n^ D '^) x 

y^McJV* /3eP(M) Ae/3 «eA 

i€N\M J 

For any /3 G V(N) with > 2 we can realize a term of the above expression, letting 
N \ M be the element of /3 containing 1, and letting 7 be the restriction of to M. 
Therefore 

K x (e 4> H(X D ,v l }) = E x ( Y f M(e )6(|/3U,6)n N ^^ri( xD '^))^)- 



I0I>2 



□ 



Lemma 2.4. Assume (Al). Let D be a domain in ~R d satisfying sup x&D E x (t£,) < 00, 
where is the exit time from D for Brownian motion. Then there exists A > such that 

supN x (exp(A(X D , 1)) - 1) < 00. 

xeD 
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Proof. Let c n = sup xeD N X ((X D , 1)™). As a\ > 0, we have c\ < 1 by (c) of Lemma [2721 
By Lemma [2.31 with <fi = we have the following recursion relation, for n > 2, 



(2.10) c n < 



E x { f ' K\P\,Q,Zt)J\c\A\dt 

/3eV(N) ^° A&p 
W>2 



(2.11) < K m \P\U c \ A \ 



/3eV{N) Ae/3 

l/?l>2 

n , j 



(2-i2) = kY'-^ y — -tk 

v ; ^ j\ ^ i 1 \i 2 \...i i \ 1 1 k 

3=2 J ii,j 2 , ...,%•>! J i=l 



h+i2+---+ij=n 



where if > 0, m<i = 2a 2 + J °° r 2 7i(dr) and = J °° r k n(dr) for fc > 3. For N > 1, let 
g-jv : [0, oo) -> [0, oo) given by 



c„ A™ 



(2.13) ^(A) = V^-, A>0 



ni 

n=l 

Using ( 12.10p . by an inductive argument it is easy to see that c n < oo for all n > 2, so g^ 
is well defined for all N > 1. Set "4>{u) = ^2^2 ~]\ u ^ = i } ( u ) ~ a i u ~ a 2 u2 for < w < A . 
Then using (I2.10p again, for A > 

N n j 

n=2 3=2 ii,i2,...,ij>l i=l 

h+h+— +ij=n 

N i \ it, 

j=2 J ' l<h,i2,...,ij<N i=l 
N 



i=2 J 

(2.14) < \ + K$(g N (\)), 

provided ^jv(A) < Ao- By the assumptions on tp, tp is infinitely different iable on [0, Ao) 
with ip" > and ^(0) = = ip'{0). So we may find Xq < Xq sufficiently small that the 
line through (x ,i/}(x )) with slope 1/K is secant to the graph of g^. That is, 

(2.15) i>{x) > ip(xo) + for x G [0, x ). 

K 

Let Ai = Xq — Kip(x ) < x < A . Since 



= ^(0) >^o)-f = -^, 



we also have Ai > 0. 
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We claim that 

(2.16) g N (X)<x VAe[0,Ai]. 

To see this, observe that g^ is continuous and strictly increasing, with gjv(O) = 0. So if 
(I2.16P fails to hold, there will be a unique A < Ai with #at(A) = xq. But by (12.141) . 

xo = fi'jv(A) < A + Ki>(g N (X)) < Ai + K${x Q ) = x , 
which is impossible. 

Now just let N — > oo in (I2.16P to complete the proof. □ 



The authors are grateful to Amram Meir, who showed them how to construct this type 
of argument. For example, see [22J. 



2.4. Potentials. We will need certain results concering potentials of specific partial dif- 
ferential equations. These potentials will be used to describe the exit measure conditioned 
to hit certain points on the boundary of D. 

Let N — {1, 2, . . . , n} be as before. Then for every non-empty subset A C N let us suppose 
we are given a solution u A > to the equation |Au = ip(u) in D. For convenience we 
also set u A = for A = 0. Define 



(2.17) v A = (-l) |A|+|B|+n+1 n B and v A = ^ (-1) 

N\ACBCN <D^BcA 

We shall assume that 

(2.18) v A > for all ^ A C N. 

The example to keep in mind is as follows: for Ti, . . . , T n C dD, let 

u A (x)=N x (X D charges (J^), 

u (x) = N X (X D charges Tj for every ieA), and 

va(x) = N X (X D charges Tj but not Tj, for every i £ A and j ^ A). 

(I2.17P holds in this case, by a simple inclusion-exclusion argument (see also Lemma 5.1 
in [29] and Section 4 in [30]). 



Lemma 2.5. ,4ssnme fl2T4"|) and (I2.18p . Then 



BCTV 0^BcA AcBcN 

AV\B+§ 



V B . 
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(b) For 

(2.19) ^-^>A = -E^7r^ E IK 

3=1 C 1 UC 2 ...UC ] =A i=l 

(c) Assume also (Al). The following then holds at any point where u A < A .' 

(2.20) iA^-^^ = -f; ( - 1)J f )(0) E (-i) |A| n^(-i) IQI 

j=2 3" CiUC2-UC j =A i=l 

where <f>(u A , v A ) = ^^^T^ * 0- 

In particular, if e rX n(dr) < 00 for every A > then f)2.20p holds without 
restriction. 



Note that all terms in the sum from (I2.19P have the same sign, whereas those from (12.201) 
vary in sign. The extra conditions in part (c) arise because of the possibility of conditional 
convergence. Part (b) describes the functions we're primarily interested in, but part (c) 
will be useful for asymptotics. 

Proof. The proof of first equality in part (a) is the same as that of (a) of Lemma 4.1 in 
[30]. As indicated in Remark 4.2 of [30], the proof of the second and third equality follow 
similarly. 

We will show part (b) first. 
Proof of fflTm : 

-Av A = Yl (-l) lAl+lBl+n+1 -Au B 

N\AcBcN 

= ^2 (-iy Al+lBl+n+i ^(u B ) 

N\AcBcN 

= Yl (-l) |A|+|B|+n+1 Uu B + a 2 (u B ) 2 + /"°°(e-™ S + ru B - 1) ir(dr)\ 

N\AcBdN \ JO J 

= olxVa + 2a 2 u N v A - a 2 ^ v c v c >+ 

CUC'=A 
C,C'^$ 



POO 

+ V / ( e ^™ S + ru B - 1) n(dr) 



N\AcBcN 



In the last line we have used the definition of Va and (b) of Lemma 4.1 of [3U] (which will 
be recognized as the current lemma in the case ip(u) = 2m 2 .) Using the expression for 
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?//(■) from (12 .7p . and that of v A we get that 

-A^4 - ip\u N )v A = - (-rv A e~ uNr + rv A ) ir(dr) - a 2 E v c v c >+ 

J ° Cl)C'=A 

poo 

V (_l)l^l+l«l+n+i / ( e -™ B + ™ B -l)7r(dr) 



N\AcB(ZN 

(2.21) 



-a 2 



POO 

V V C V C , + V ^AMBl+n+l / {e -ru* + _ 1} ^ 



C*UC*'=A N\AcBcN 



Consider the last term of this expression. A ^ so by Lemma 12.11 

POO 



N\A<zBcN 

Therefore 



POO 

J- (-l)\A\+\B\+n+l / (e -™* + ru B e -ru» _ ^ 

A , — T~t , — »T JO 



N\AcBcN 



/_ 1 yA|+|B|+n+l / 

Z^t^AT JO 



{e~ ru ~ +ru B e~ TU " -ru N e~ ru " - e~ ru " ) 7r(dr) 



N\AcBcN 



poo 



N\AcB(ZN 



= (_i)W+|B|+»+i / e -™^I_(^_ n £ )%(dr) _ 

By part (a), u N — u B = E] t>c- So by monotone convergence we have that the above 

cns=0 

is 

N\AcBcN j=2 J° J' %^CdN 

cns=0 

OO /*00 7' 

(2.22) = W ^3 E (-l) |A|+|B|+ " +1 ( E ^c)^(dr). 

3=2 ^° N\AgBgN %^CdN 

b^% cns=0 
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Using standard multinomial expansions and Lemma 12.11 we observe that j-th summand 
is, 



^2 (-l)\ A \+\ B \+ n+1 ^ ^ ■ 

N\AcBcN 9^CcN 

CDB=<D 



J 



N\AcBcN C u C 2 ,...Cj i=l 

%^C l <ZN\B 

= E (]><*) E 

d,C 2 ,...Cj i=l N\AcBcN\\J.Ci 

= e (n^)(-i) |Ai+ra+i (-i) n " iA| i A ^ Ci 

Ci,C2,---Cj i=l 

= - e rW 

CiUC 2 ...UC* J = J 4 i=l 

Together with (I2.2ip and (I2.22p this implies (I2.19p . It is clear from the proof that the 
sum in ( 12 . 1 9[) converges. 

Proof of (1^201) : 

l -Av A = E (-1) |B|+1 ^A^ 

%^B(ZA 



E (-i) |B|+ V( 

%+B<zA 



Observe that u B < u A by part (a) and (IH8| . So by (|2~8|) . we can expand tp in a series, 
at any point where m a < Aq. Therefore 



(2.23) = aX + E ( ~ 1)3 'f )(0) E 
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For j > 2, we have by part (a) that 



0^BcA 



E E n^c- 1 ) 1 * 1 

0^BcA fr£Ci,C a ,...Cj<zB i=l 

£ n^(-l) |Cl| l(A = Ui 1 C7 l ) 

<t^C 1 ,C 3l ...,CjCA i=l 

(using Lemma 12.11 ) 

(-i)i^ e n^c-i) 1 * 

i=l 



0^Ci,C 2 ,...,C^A 



fe=l ^ ' \0^Ci,C 2 ,...,C,_ fc ^i=l 



4) W+i ^ J^(-l)l^l 



i=l 



A=ul =1 d 



A=\J i=1 d i=1 
$£C\,Ci,...,C i ±A 



J 

A\j 



E E ) +(-i) |A|+1 (-w A ) j -(-w A + -(-i) |A 

0^BcA \%jLC(lB,C^B 



V 



(2.24) 
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Using ( 12T23|) and f[2T2ij) we have 

Iaz/ = 0l ^ + ( -i)^-»g; tWp® yj n^(-i) |c '' + 

f; ( ~ 1)J f )(0) (-n^ - (-^ + -(-i)^i^y 

i=2 J ' 



= a x+(-i)' A ' +i f: ^^ e n^(-i) iai + 

<b^Ci,C2,-,Cj^A 

POO 

JO 

0^Ci,C 2i ...,C^A 



which gives (I2.20p . Positivity of <f> follows from monotonicity of if)'. □ 



3. The Branching particle description 

In this section we shall define a martingale change of measure which will represent various 
conditionings of the exit measure of if} super-Brownian motion. For each such condition- 
ing we shall also present a branching backbone representation, in which the conditioned 
process is realized as an unconditioned superprocess with immigration of mass along a 
branching tree. 

Let ^ A C N = {1, 2 . . . , n} and let u A > and v A > be as in Section E2D Specific 
examples will be described later, but even at this level of generality we can use these 
functions to define an associated martingale. Define 

oo m 

(3.1) M fc = exp(-<X fc ,^>)]T- n<**.«<*>>0. 

m=l ' C 1 U...UC m =N i=l 
C#0 



In Lemma 4.3 of [30], it is shown that this sum is finite, and 

M fe = £(-l)Wexp(-(X* 

AcN 
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(note that the A = term equals 1, by our convention that vr = 0). 

From this or otherwise it can be easily checked that is a martingale under N x , with re- 
spect to the filtration T\.. So for $ measurable with respect to F k let M x . = v J-^ N x (&Mk) 
be the associated Girsanov transformation - that is, the /i-transform, for the "harmonic" 
function 

h(jM)= ^(-l^'expH/iV)). 

AcN 

Dynkin has developed a general framework for such /i-transforms of superprocesses, which 
he calls X -harmonic functions. See [IT] . 

Branching backbone : We describe the direct construction of M x in terms of a 
branching particle system tracing out a backbone along which mass gets created. Let 
Lf = L^i ou n f = tjAJ — ip'{u N )f. Essentially we will constructively generate a measure 

N x and then show that it agrees with Wl x . N x will have a branching backbone T equipped 
with mass creation/immigration both along the branches and at the nodes. Each branch 
will be labeled by a nonempty subset A of N. The mass created along the backbone 
will evolve as an unconditioned superprocess but with a modified killing rate and Levy 
measure. The first order of business is to give a more precise description of the various 
components of this process. 

Evolution of Mass: Once mass has been created, it evolves as an unconditioned super- 
Brownian motion but with a modified branching law. More formally, for $ measurable 
with respect to let N y ($) = N 2/ ($e^ JV ). We will show that this measure indeed 
describes a superprocess. 

For any A > and y G D define $(y, A) = ip{u N {y) + A) — ip{u N {y)). It is easily checked 
that 

poo 

tf( y , A) = a x (y)\ + a 2 \ 2 + / e - uN{y \e~ Xr - 1 + Ar) 7r(dr) 

Jo 

where di{y) — a± + 2a 2 u N (y) + J °°r(l — e~ ruN ^) n{dr). This is therefore a spatially 
varying branching law, of the form (12. 3p . But now in places where u N is large, we will 
have that the killing rate d\ becomes large too (provided a 2 ^ 0), and the Levy branching 
measure 7i(dr)e~ ru becomes small. 

The following is a special case of Dawson's Girsanov formula (see Theorem 7.2.2. in [7]), 

Lemma 3.1. Assume (12.41) . Let N x be the excursion law for super- Brownian motion in 
D with branching mechanism ip. Then for every <f> > 0, 



N x (l-e^ xk ^)=N x (l-e-( xk ^). 
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Proof. Let g be the solution to 

^&g(y) = 4i{y,g{y)) for y e D k 
g = (f) on dDk- 

Then f = g + u N is the solution to 

±A/ = r/>(f) in D k 
f = (f) + u N on <9-Dfc, 

so 

^(1 - e"^) = N*((l - e -(^^)e-^^)) 
= N z (l - e -<* fc .^>) - N.a - e^ xk ^) 
= (g{x) + u N {xj) - = <?(x) = N x (l ~ exp -(X fc , 0)). 

This suffices. □ 

T Backbone: Under N x we begin with one particle in the system that performs a vn 
transform of the motion with generator L in D. Note that by Lemma 12.51 is L- 
superharmonic (as indeed are all the va)- Let — Va(-) denote the right hand side of (b) of 
Lemma 12.51 

If the particle dies at some site y, a random number of particles are born in its place, and 
numbered 1, 2, . . . , j. The first particle is then assigned a randomly generated tag A±, 
the second is tagged with A2, etc. For j > 2, the probability that j > 2 particles are born 
and that their tags are a specified sequence A±, . . . , Aj is given by: 

(3-2) 1 , b(j,u N (y),y)f[v Ai (y). 

Here the Ai are chosen so that Ai 7^ and U{ =l Ai = N. In full generality, particles need 
not die before exiting D, but in our main example below they will in fact always die in the 
interior of D (and the va will accordingly be //-potentials). Note that the above defines 
a probability measure, by definition of Vn- 

Each of the particles so created now evolves as a va transform of the motion with generator 
L (where A is the particle's label), until it dies. Whereupon a new collection of random 
branches is created, each labeled by an A\ 7^ with UjA'j = A, etc. 

It will be convenient to describe this backbone as follows. Let n t denote the number of 
particles alive at time t. Label them using an index i, 1 < i < n t , and for each one let 
Xi(s), < s < t denote the location of the particle or its ancestor at time s. Let be 
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the measure putting unit mass at each point Xi(t) for which Xi has not exited by time 
t. Then T k represents the backbone killed upon leaving D%., and we recover the whole 
backbone by letting k — > oo. Without comment we will feel free to refer to T fc in terms 
of the underlying particles, though formally it is still a measure- valued process. 

Immigration at nodes: The birth of j particles at y is accompanied by site-dependent 
creation of mass. Conditional on the location being y and the number of particles being 
j > 2, the mass created is a random variable Rj > whose law fij y (dr) is given for r > 
by 



Hi, y {dr) 



r j e -ru N (y) w ( dr j . 

2a 2 +/ ° o r 2 e-'-" JV W7r(dr)' 3 ~ ' 



In the case j = 2, also has an atom at r = of size 2a2/[2a2 + / °° r 2 e~ ruN ^ ir(dr)]. 
In other words, if a2 > then it is possible for no mass to be created at a branch with 
j = 2. In any case, all mass created at the node then evolves according to N y . 

Temporarily writing Y^ y for the exit measure from resulting from such a creation of 
mass, and Ej tV for expectations under this conditional law, we can write Y* y as J v N(dv) 
where N is a Poisson random measure with intensity n[dv) = Rj'N y (X k G dv). Thus 

(3.3) = E hy [e- R ^ l -*- {xk ^} = [ e ~rN y (i-e-(- k ^) Hy{dr) 

J[0,oo) 

b(j,u N (y)+N y (l-e^ xk ^),y) 
b(j,u N (y),y) 

Denote this quantity by M k (j, (f), y). 

Immigration along branches: 

For any A > and y € D define 

V ( y ,\)=ip'(u N ( y ) + \)-i/(u N (y)) 

= 2a 2 \+ [ re- ruN{v) (l-e- Xr )ir(dr). 
Jo 

Notice that this has the same form as ip' in (12. 7p but with a\ = and a spatially varying 
n y (dr) = e- ruN ( y) n(dr). 



We create mass along the branches. The mass created to time t forms a spatially depen- 
dent Levy process, with Levy exponent rj. In other words, if L t is the mass created until 
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time t, then 

n x (e~ XLt ) = N x {e~^^ iy ' x)T " {dy)ds ). 

This mass then evolves according to N y . In particular, if r denotes the first branch time 
and £t the position along the first branch to time t then 

If Yq t denotes the exit measure of the mass created along this branch till time r then a 
calculation similar to (13. 3p shows that 

(3.4) N x (e-< y B^>) = ^ (e -/;%(i-e-^ fc .«)dL t) 
In other words, we obtain an expression similar to that of (12. 5p . 

With this, we've finished describing the construction of an exit measure Y k under a 
probability measure N x , and are ready for the following: 

Theorem 3.1. Assume conditions (12. 4p and (I2.18p . Then 

M z (exp-(X fc ,0» = N z (exp-(r fc ,0>). 

Note that this shows that, under the two measures M x and N x , the exit measures from 
each Dk have distributions which agree. Using historical processes (as in [29]) one can 
show that the same result carries over to the process level. That is, that the laws of the 
full superprocesses agree under these two measures. 

Proof. In the present context, it is useful to label all the particles of T k that exit D^, by 
placing an order on them. So let Fk be the set of such particles, and set 7fc = \Fk\. For 
Ac N, let 

S m {A) = {(d, . . . , C m ) : C x U ■ • • U C m = A, ^ Q Vz}. 

If 7^ = m, choose at random an ordering of F^, and for A = (Ci, . . . , C m ) G S m (N), write 
T fc « A for the event that the ith particle is tagged with the set Ci, % = 1, . . . , m. Thus 
for example, 

(3.5) n x ( lk = m)= ^(T fc ^A). 

Note that if M\, . . . , Mj are disjoint ordered sets, with \Mi\ = hi and m = Yl^i then there 

are kA orderings of S = IjMj which are compatible with the given orders on each 

Mj. In other words, if a is any order on S, and if £ is an order on S picked at random, 

then the conditional probability 

. . fci !fc 2 ! . . . kA 

3.6 P(E = er E Mi = a Mi ,i = 1, . . . , j )= 3 

ml 

(writing o M etc. . . for the restriction of a to M). 
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As described initially, the root particle of the tree is always a v ^-particle. It is convenient, 
for purposes of induction, to allow the same notation to cover the situation that we start 
with our root being a w^-particle for some A C N. In this case, (13. 5 p still holds, but with 
A G S m (N) replaced by A e S m (A). With this in mind, we may define another restriction 
operation as follows. For 1 < z'i < • • • < < m, set 

(Ci, . . . , C m )|{j 1; ... )ifc } = (C i]L , . . . , CjJ. 

Thus, if A = (Ci. . . . , C m ) E S m (A) and M C {1, . . . , m}, we will have that A| M e S m (B), 
for i? = UjgjuCi. As a shorthand for the latter, we write A(M) = UigMCj- 

The case of interest is that the first branch of T fc partitions {1, . . . , m} via the descent 
relation. If (3 is this partition then there are \/3\ particles born at this branch, and 

(3.7) \j3\l ways of tagging the 1st particle, 2nd particle, etc. 
with distinct elements of /3. 

We will show, by induction on m > 1, that for ^ A C N, and (C±, . . . , C m ) G S m (A), 

(3.8) N z (exp -(F fc , 0) >T k « (d, . . . , C m )) 

= ^^T^(eJ +ttJv n( Xfc '^))- 

Taking A = A r and summing over S m (N) will then establish the theorem. 

We start with the case m = 1. Here T fc will have a single wat- process with lifetime £ > r^. 
Therefore, 

ft a .(exp-<y*,0);T* « AT) 



rpl/j' ou N ,V n 
X 

1 



x C<r fc e 



Vat(xJ 

1 
1 

Vat(x) 

1 



-.EL 



E x [v N (£ Tk )Af Th (e 

<l>+u N )] 

R x (el +uN (X k ,v N )), 



where the second last equality follows from definition of r\ and J\f, and the last equality 
follows from Lemma [2.21 (c). 



Turning to the inductive step, let m > 1 and assume the inductive hypothesis for all 
A G N, and for all values smaller than m. For simplicity, we will verify (13. 8p in the case 
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A — N. For ( the lifetime of the initial particle, and A = (C±, . . . , C m ), we have by ( 13. 4p . 
(E3D, (EZD, (E2D, (ESD, and the definition of T, that 



N^e-^'^T* « A) 

= E # 



l C<7 >e -J ° »»(e».N«.(l-eJ))cfa x 



p€V({l,...m}) 

|j8|>2 



xM(|/3|,0,£ c ) x x VUgl! X 

x n ^, |A|! x 17^(6"^; T*«A| A ) 



By (T22D and (JSHJ) this equals 



m!u7v(x 



E 



0GP({l,...m}) 
V l/3|>2 



-/o'7(6.N €s (l-efc))d S>< 



1 



xb(\fi\, (f> + u N ,t t ) x J] |A|!t; A{A) (6)%(e-^^;T fc « A| A ) ) eft 



/ 



-.EL 



/3eP({l,...m}) 
\ l/3|>2 



■J?[^ / (tt JV (6))-Hjtt.^ C ,(l-e*))]d. x 



X 



b(\(3\,cj> + u N ,t t ) J] |A|b A(A) (6)%(e-^);T fe » A| A ) ) dt. 
Ae/3 



By definition of r\ this is 



■E x 



E r M(e^)6(|/3U + ^,6) x 



/3£P({l,...m}) 
\ l/3|>2 



x J] |A|!t, A(A) (6)%(e-< yfc ^;T fc ~ A| A ) J dt, 

AefS 
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-En- 



mlVN(x) 



^2 r ^H + ^)K\p\A+u N ^ t )x 



/3eP({l,...m}) 
\ l/3!>2 



Ae/3 



H(x k ,v Cl ))) dt. 



By Lemma [2.3I we conclude that 



rn 



i=l 



□ 



4. Conditioning the exit measure to hit n points 



In this section we shall consider the exit measure when it is conditioned to give positive 
mass to n small balls on the boundary of D. We shall study the limit as the radius of 
these small balls tends to 0. In the first part, we set notation. In the second part, we 
discuss the case when ip is analytic. Here the limit is given by a martingale change of 
measure as in the previous section. Then we describe an "explosion" phenomenon when 
n = 2 and = A 1+/3 , < f3 < 1. Finally we discuss instablity in the limiting process by 
establishing a range of possible limits. We begin by fixing some notation for this section. 



4.1. Conditioning to hit n small balls. Let D be a C 2 -bounded domain in IR d , for 
d > 4, N = {1,2, ... ,n} for n G N, and let {zi, Zi, . . . , z n } be distinct points on the 
boundary of D. Let U be the potential operator for Brownian motion killed upon hitting 
the boundary of D. For x G D let K®(-, •) be the Martin Kernel on D and let G £>(■,■) be 
the Green function. For any z G dD and e > we define A(z, e) = B(z, e) fl 3D where 
B(z, e) is the Euclidean ball of radius e in R d . 

For any A G N, set 

B? = U i€A Bl = U ieA B{ Zi ,e) and = U ieA A\ = U ieA A(z h e). 
Fix < 5 < 1 such that B(z iy S ) fl B(zj, S ) = if i ^ j G N and for any e > 0, set 



= D\U ieA B( Zi ,e) and r e = = t d a. 
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For any x G D, define the functions: 

vt{x) = N x .(n(^ D ,lA(^))>0), 

i€A 

and 

v a,e{%) — N s ( X D charges all A* for i G A and does not charge for j G" A ). 
It is easy to see by an exclusion-inclusion argument that 

N\ACBCN 

< = - E (-^'V, 
(4.i) ^ = - E (-i)'V- 

Thus f[2~TTj) and f[2~T8]) hold. By (a) of Lemma [2J3 |Aw^ = V(^). 
We will need the following two Lemmas in section 14.21 

Lemma 4.1. Let D be a bounded C 2 domain in dimension d > 4. Let B be d-dimensional 
Brownian motion started from x G D, under a probability measure P x . For y G D let rriy 
be the harmonic measure starting from w. Let z G 3D and z G 3D. Then for x G D 
fixed, 3 constants Ck, k = 1, 2, 3, 4 st/c/i £/ia£ 

(a) P y (B T2e G 0B(*,2e)) < ^(Af^)) for y G D \ B(z, 4e), 

(b) f D i/XCA^, e)) 2 rfy < L7 2 e 2 mf (A(z, e)), 

(c) (y, z) <C 3 \y- z\~ d dist(y, 3D), and 

(d) m x (A(z,e) >C A e d ~\ 

Proof. This uses the proof of Lemma 5.4 in [29], which in turn uses an argument from pp. 
More specifically, (a) follows from (5.6) in [29] and (b) follows from (5.8) in [29]. We note 
that in [29] the domain D is assumed to be bounded Lipschitz domain which includes the 
class of C 2 -bounded domains. Part (c) follows from Lemma 2.1 in [TJ. Part (d) is just the 
fact that the density of harmonic measure with respect to surface area is bounded away 
from 0. □ 



Set p = (d — 3)/(d — 1), and note that < p < 1 since d > 4. 
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Lemma 4.2. Let D be a bounded C 2 domain in dimension d > 4, A C N = {1, . . . ,n}, 
and A > 0. Then there exists 9 > and e > 0, such that whenever < e < e and 
y G Dg eP , V % G A then 

v\{y) < A. 



Proof, (a) By the maximum principle, it is enough to prove the lemma for y G Df eP and 
dist(y,aD) < Be?. 

Fix an i G A, and denote v\ by v t and Zi by z. Let 6 1 > 2, e < min(5 , 1),0 < e < e . 
Using the comparison principle (see 8.2.H in jTUJ ) , and (2.2.10) in [32] there exists C\ > 
such that 

(4.2) v e (x) < ci dist(x,9D)- 2 ,Vx G L>. 
Now using the Feynman-Kac formula we have, 

Ve(y) = E y L(B(T2e))eM- ^f^f ds ) 1 ^ < ^)) 

< S„(i; e ( J B(r2 e ));l(r2 e <r D )) 

< Cie~ 2 P^(r 2(; < t d ) 

(4.3) < c 2 e- 2 P y (B TD G A(z,e)), 

where the last inequality follows from the proof of Theorem 3.1 in [Tj. As p < 1 and 
9 > 2, < e < ^\y — z\. Using the above and Lemma 2.1 in [I], we have that 

(4.4) v e (y) <c 3 \y-z\- d dist(y,dD)e d - 3 . 
Since D is bounded we have for all y G -D#*p, 

v e (y) < c 3 (9e p y d+1 e d ~ 3 < C3 0-d+i e d-s+(-d + x) P = ^-d+i 

As d > 4, we can choose 9 large enough to obtain the required upper bound. The argument 
is independent of % and the proof of the lemma is complete. □ 



4.2. Asymptotics for analytic ip. In this section we shall assume that ip is a real 
analytic function. We also assume that ai = (a version of the results should hold 
for a\ > as well, but for the Martin kernel of killed Brownian motion rather than of 
Brownian motion itself). Note that if a± = then (Al) =>- (I2.4p . The results presented 
here and their proofs largely mirror those presented in [29] which considered the case 
= 2m 2 . The principal difference lies in the proof of Lemma 14. 3[ where a more 
delicate argument is required in order to obtain convergence of the power series arising 
there. 

Theorem 4.1. Let D be a bounded C 2 domain in dimension d > 4. Assume (Al) and 
thatai = 0. Let A C N = {1, . . . , n}. Then for x,y G D , 

lim vf(y) k a {) 
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where 

K?(y,z t ), A = {t}, 



As part of the argument, it will emerge that K A is actually finite. Note that if \A\ > 2 
we may also write 

(4.5) tfX^gLM?) E U^K^y). 

Ciy^dt and disjoint 

It is the latter form that comes naturally out of the formulae of Sections [2] and |3j 
We need several lemmas before starting the proof. 



Lemma 4.3. Assume the conditions of Theorem \4-l\ and let p, €q, and 9 be as in Lemma 
4-2\ with A < Ao/2 n . Then there is a C < oo such that VO < e < eo and y G Df^OeP 



Proof. We prove the lemma by induction on the size of A. 

Step 1: = 1) Let A = {z} and v e = vf. As in the proof of Lemma \4.2\ we use (14. 2p . 
the Feynman Kac formula, and (a) of Lemma 14. 1[ to get 

ve(y) = E y (v e (B(r 2e ))eM- jT ds)j 

< E y (v e (B(r 2e )); B(r 2e ) G dB(z, 2e)) 

< Cl e- 2 P y {B(T 2e ) G dB{z,2e)) 

(4.6) < c 2 e- 2 mf(A( 2 ,e)) 

for y G Df t . Using the Palm formula (Lemma 12.31) with = and n = 2 we then obtain 
that 

(4.7) N x (X D (A(z,e)) 2 ) = ^ 2 \0)E X / [N Bs (X D (A(z,e)))] 2 ds. 
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Using the above, the Cauchy-Schwartz inequality, Lemma 12.31 and Lemma 14.11 (b) we 
have 

v e (x) = N x (X D (A( Zl e))>0) 
N x (X D (A(z,e))f 
~ N x (X D (A(z, e)) 2 ) 

mP (A(z,e)) 2 
^)(0) J D G D X (x,y)m°(A(z,e)y)dy 

(4.8) > c 3 e" 2 mf(A(z,e)). 

Consequently, (14. 6ft . (14. 8ft . and the boundary Harnack principle (see [3]) yield 

v e (y) e*m°(A( Zl e)) 

for y G D± e . By decreasing e if necessary, we can (and will) assume that Di d D^ tP . 
This establishes the case \A\ = 1. 



Step 2: (\A\ > 1) Assume the result for every proper subset of A. Set q — \A\ and 
a(k) = 2 k 9e p for k = 1, . . . , q. By hypothesis, each v\(y) < A /2 n on D^, X y The same is 

then true of v B for B C A, so by Lemma [2.51 u A < Ao on Thus (12.201) applies and 

by the Feynman-Kac formula, we obtain the following on D\ : 



EJe- 



{Br)dr vt{B T \) 



J'=2 



*(«-!) 



Juf'(S t )(-l)l c< le-/o^)*dt 



i=i 



(4.9) 



/ + //. 



Here < 4>f = (f>(vy, uf). Consider the first term in 04. 9p : 



I < E y (vt(B Taiq _ 1} )) 

< Evi^Br^HB^ edB( Zi ,a(q-l)))) 

< SU P v*(-) PyiB^EdB^aiq-l))) 

< c 6 J2 SU P vA(.)m°(A( Zi ,a{q-2))), 
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where the last inequality follows from Lemma [4.11 (a) [with for z and a(q — 2) for e]. 
We know that vf < vf^ l \ Using this and (14.81) with z = Zi, we have 

I < v-^ suPaB( 2< , a(q -i)) vf" {l} {-) m°(A( Zi , a(q - 2))) 

u , \- sup aB ( z . ia(g _ 1)) u e AU>} (-) m^(A(^, q(g - 2))) 

1 ' j - C7 ^ rw^(*) - 2 <(A(,, e )) • 

By the boundary Harnack principle, 

m°(A{z h a(q - 2))) < c 8 mf (Afc, a{q - 2)))K x (y, Zi ), 

and as in mf (Afo, a(g-2))) < c 9 a(g-2) d ~ 1 < c 10 e d ~ 3 . Likewise e~ 2 mf( A e)) > 
cne d_3 , by (d) of Lemma 14. II By induction, '(■) / 'FLeA j^i v e( x ) is bounded, in partic- 
ular, on D^^s by cu Y2jeA\{i} ^("' which by Lemma [47T1 (c) is bounded by c\^a{q — 1) 
on dB(z,i, a(q — 1)). Therefore 

(4.11) - 1 <c 14 «(g-l)^^(y,^) 

Let 

3 

E j = {{(';}> . : ^ Q ^ A U^ =1 C, = A, and \A\ + \ C i\ is even ) 

i=l 

and 

j 

O j = {{ay i=1 -.Q^Ci^A, uJ =1 Ci = A, and | A| + ^ |a| is odd}. 
Then the second term (II) in (14.91) is, 



i=l 



j=2 J - - /n i=l 

-i'-"T""E^ pn^-" 

j=2 ^' Oj i=l 

j=2 "'' E,- 1/0 i=l 

as all terms in the second summand are non-negative. Therefore 

(4.12) // < f- Ma^p^ 
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First consider the case j < \A\ and {Ci} J i=1 G Ej. We observe that for k ^ I, if Ck\Ci ^ 
then 

yCkyCl v C k \ClyCl 

So every term in this sum with j < \A\ is dominated by another term in which {Ci]\ =l G 
Ej are such that C k PI C\ = for k ^ I. Thus it suffices to bound such terms. For disjoint 
{Ci}U G E h 

/•ra(,-i) rri r«(.-i) FH . Y\ r K x (B t ,z k ) U,ar vUx) 

lL e A^) " Jo UieA^x) 

(4.13) = 0E y j ^ ^ i^(5 t , * fc ) di. 

^° i=i keCi 

Let fci, . . . , kj be distinct. Then 



E y ]J i=l K x (B t ,z ki )dt<G D (f[K x (;z kt ))(y) 

Jo i=i 

j j 
< c 15 ^2G D (K x (-,z ki ))(y) <c 15 ^2K x (y,z ki ) 

i=l i=l 

\A\ 

(4.14) < c 16 J2 K *(y,zk) 

k=l 

where the second inequality is due to the fact that z^s are separated by So and the third 
inequality follows from the "3-G" theorem (see (5.17) [29]). It follows from (14.131) that 

(4.15) Ey n 1 :L dt ^ c ^ K ^^- 

J0 LLiGA V e\. X ) k=1 



Now consider the case j > \A\. For {Cj}f =1 G Ej we can select \A\ of the whose union 
is A, and apply the above bound to them. For the other Ci we apply Lemma I4.2[ which 
gives that v^ 1 < m&x k < A, where 2' A U < \ . Using ( I4.12p . and (I4.15P we have 



II 



< Ci 7 



< C17 



^i-iy^(o)\Ej 

j=2 



i=2 



j=|A|+l 
oo 

- £ 

j=|A|+l 



-1) J >^(0)|^|A J - 
j! 



■14 



14 



^2K x (y,z k ) 



k=l 



(_l)^(i)( )2Jl^lA^I A l 
j! 



14 

^2K x (y,z k ) 

k=l 



S Cl8 z k ). 



k=l 



□ 



30 SIVA R. ATHREYA AND THOMAS S. SALISBURY 

Lemma 4.4. Assume the conditions of Theorem \4-l\ Let y e Df, where 5 < S . Then 



uniformly in z G dD, 

limP^(exp(- [ S <f)f (B t ) dt)) = 1. 
£ ->° ' Jo 

Proof. We have v\{-) — > as e — > 0. So by the assumptions on and Lemma [4.31 and 



-1)1-11 



Consequently it suffices to prove that 



-»■ 0. 



limPj /z (exp(-Ar 5 )) = 1, 

A^O 



uniformly in 2. As .D is Lipschitz, we have 

sup P yz (r s ) < oo 

yeD s ,zedD 

by Then Lemma 3.7 in [I] implies the result. □ 

Proof of Theorem \4-l\ We will use induction on the size of A. First consider the case 
where our target is a single point. Let A = {i}. When convenient in the proof we will use 
Ds for Df for some 5 > 0. Let x, y G D<5 . By the Feyman-Kac Formula, for each fixed 
5 < 5q we have 



vj(y) E y (v\(B Ts )eMr 5 <Pe(Br)dr)) 
vi(x) E x (vl(B T6 ) exp( J T * #(B P ) dr) 

/ flgf P„(exp(- g #*(g t ) <ft))lff '(y, z)^(z) m*(dz) 

f dDs P xz (eM~ J? <tt(B t ) dt)X(z) m^(dz) 
J dDg P yz (exp(- J TS tf(B t ) dt))K?'(y, z)v\{z) m°*(dz) 
j Ds vi(z)m x 5 (dz) 
J Dx vl(z)m^(dz) 

x 



The measure 



J 9Dg P xz (exp(- Io S WW dt))vi(z) wS'{dz) 
v\{z) m x s (dz) 



f vi(z)m x s (dz 



is a probability measure on dDf. Since the boundary of Df is compact, by Prohorov's 
theorem any sequence €j has a subsequence, again written ej, for which X ej> §(x, dz) =>■ 
Xg(x,dz) weakly in the space of probability measures. Also, K x s (y,z) is continuous and 
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bounded in z, for z G D fl dDf, when x, y G Df Q . Consequently, Lemma [4.41 implies for 
ij6 Z?^ and for all 5 < 5 



K\y) r 

lim ^TT = / K^(y,z)X s (x,dz). 



'dD 5 

The limiting function is harmonic in y for y G Dg . By a diagonlization argument, we can 
assume there exists a convergent subsequence of our sequence such that the convergence 
holds simultaneously for a sequence of 5/s which converge to 0. By Lemma 14.31 we see 
then that the limit is harmonic in y with boundary value on 3D n dDf for all 5 > 0, 
and is 1 at y = x. This implies that the limit is the Martin Kernel for Brownian motion 
in D. This all subsequences have a subsequence which converges to the Martin kernel, 
and so the limit itself exists. 

To prove the induction step, fix A and assume that the result is true for all proper subsets 
of A. Therefore if Uf =1 Cj = A and ^ d ^ A we have 

FE=1 V e l {v) [ H*=l K ^(2/) 1 {^nC i =0, X<k#<n} ^ 1 < 3 < \ A \ 

riieA^^) o otherwise. 



Let x, y G D. Let e > 0, 77 = 9e p > be small enough so that x, y G -D^ and Lemma I4.3[ 
Lemma |4.4[ and Lemma [4.21 apply. By the Feynmann Kac formula and ( I2.20p . we have 



then 



^) > (-1)^(0) V£ ri^We-j?^)*^ 



By Lemma 14.31 Lemma 14.41 the dominated convergence theorem (which applies as in 
Lemma 14.31 by the bound provided by Lemma 14. 2[) and the induction hypothesis we have 

,4,6) ta wJ% > £<ZM E B Pf^)*. 
For the upper bound, as in the proof of Lemma 14.31 we have 

vf(y) < E(v?(B Tn )) ^ (-1)^(0) p nti^W ^ 
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Letting e — > 0, the first term — > (since r) — > 0). Using the inductive hypothesis and 
dominated convergence theorem for the second term we obtain 

«-0 ll i6 A» 6 W i=3 J- U 3 i _ 1 d=A J ° *=1 

Ci^0=C fc nC; 

(l43j) . ( ESP and f HTTj) now yield the result. □ 



Theorem 4.2. Assume the conditions of Theorem \4-l\ Let G J-fc be bounded, and fix 
x G D. Then for y G D, 



JimN y (^| J]^, 1 A .) > 0) = J ±—K v (* k MS), 



where 

<= e n 

Proof. We will need two preliminary Lemmas. 

Lemma 4.5. Assume the conditions of Theorem \4- 1 . S'et W 7 ^ = exp((— l^^X*, )) — 1. 
TTien 

N,fn(^,i A ,)> |^) = £ (n 1^1) (-1)^-^' 

\t=l / j=l CiUCi-Cj \i=l J 

CiUC 2 -UC 3 =N 

Proof. Using the arguments presented in [9] or otherwise one can verify that 

<^) = N,(l-exp(-A(X B ,]Tl A »>)) 

increases to uf(y) for all y G D as A — > 00. From here on, the proof of this lemma is the 
same as the proof of Lemma 5.8 in [22]. We will not present it again here, other than to 
remark that it is based on the Markov property of exit measures. Note that a different 
indexing system is used in [29] , which accounts for the j\ factor. □ 



Lemma 4.6. Assume the conditions of Theorem \4-l\ Let G be bounded. Let 
Ci, . . . , Cj be distinct and nonempty, with U J i=1 Ci = N. Then 

So FT vHx) = ^(*fcll\ X i K *)) l { Cl disjoint}- 



Proof. The proof of this lemma can be obtained by imitating the proof of Lemma 5.9 in 
|29j . The only changes one has to make are to use: Theorem 14.11 in place of Theorem 5.3; 
Lemma [4.31 instead of Lemma 5.4; Lemma [2.41 instead of Lemma 2.7. □ 
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To complete the proof of the Theorem, observe that 



N„($ fc | l[(X D ,l Al ) >0) 



8=1 



N,(nr=i(^iA l )>o) 

^(^N,(nr = i(^ D ,iA>) > o|j*)) ^ niu^) 



X 



By Theorem 14.11 



By Lemma 14.51 and Lemma 14.61 



^1 ^ N tf ($*nti IW^K-lf+EL^h 



i =1 j! c.ufe.c, IE=i*) 

CiUC 2 ...UC* j =Af 



2^-1 

— > 



E -i E N »(** II< Xfc > disjoint } ) 



j=l J dUC 2 ...Cj i=l 

C 1 UC 2 -UC j =N 



= e n y n = ^(^O- 

crG-P(A) CGo- 

The statement of the Theorem then follows easily. □ 



4.2.1. Weakening Hypothesis (Al). It should be possible to weaken the assumption (Al) 
in Theorem 14.11 to the following: 

(A2) f™r n ir(dr)<oo, 

when TV = {1, 2, . . . , n}. We are able to prove Theorem 14.11 under this weakened assump- 
tion, though only in the case n = 2. We begin with a lemma. 

Lemma 4.7. Let D be a bounded C 2 domain in dimension d > 4. Assume (A2) and that 
ai = 0. Let e > 0. There exists a, b G E, C\ > 0, f e , g e , h e : D\ )2 — » [a, 6] swc/i t/iat /or a// 
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y e D]>\ 

(4.18) v\\y) > E y (e~fo a M*>*v?{B Ta )) + E y ( f" g t {B t )vl{B t )v 2 e {B t )) dt 

Jo 

(4.19) v^(y) < ^(e-JiT , /.(*)V a (B T J) + ^(r\(^K 1 (5t)t; e a (B t ))(ft ) 

Jo 

(4.20) ]img e (y) = lim/i £ (y) = ip {2) (0) 

Proof of Lemma: Observe that, 
\av^ = £ (-l)'*' +1 ^f 

0^BC{1,2} 

= E (-i) |fl|+ V(«f) 

0^BC{1,2} 

= ^) + ^(v 2 ) - ij)(v\ + v 2 - v 1 /) 

= -[i>(vl + v 2 ) - + Mv 2 e ) - V(0)] + VP{v\ + vl) - 4,{v] + v 2 - v 1 /)] 

When n — 1, (A2) implies that ^ is twice continuously differentiable. Using Taylor's 
theorem on ip, we have 

Uv 1 / = -[^(vl)v 2 + ^M^f] + [^(0)v 2 + ^ 2 \a 2 ) { 4f} + [^("sK 1 ' 2 ] 



2\2 



-i^\vl) - ^ {1) (0)k 2 - l¥ 2) M - ^ 2 \a 2 )] y -^- + ^("sk 1 ' 2 



2\2 



-^ 2 \a,)v\v 2 e - [^)( ai )-^ 2 \a 2 )} K -^ + ^ 1 \a 3 )v 



1,2 



where ojj : .D*' 2 — >■ [0, oo) are measurable functions such that v\ < a,\ < v 2 + v},0 < a 2 < 
v 2 , v] + v 2 - vl> 2 < a 3 < v 2 + v\, and < 0:4 < v\. 

Repeating the above calculation with the roles of vl and v 2 reversed, we would obtain 

\avI> 2 = -^( 2 >(/5 4 )^> e 2 - [^ (2) (/30 - ^ (2) (/5 2 )]^ + ^ (1) («s)^' 2 , 

where fa : D^ 2 — > [0, 00) are measurable functions such that v 2 < Pi < v 2 + v\, < fa < 
vl, and < ^ 4 < v 2 - 

Using the Feynman-Kac formula, we have 

vl> 2 (y) = E y (e-{° a ^ W ^ B ^ dr vl> 2 (B Ta )) + E y ( H e~ ^ W ^ B ^ dr J t {B t )) dt, 
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where J t : D — > R is given by 

J £ = V (2) KKV+^( 2) (ai)-V (2) («2)) 



(v 2 ) 2 



2 

1\2 



= ^ (2) (/34)^ 2 +(^ 2) (A)-^ (2) (/3 2 ))%^- 
Now, ^ 2 (-) is a non-negative decreasing continuous function. If v\ < v\ then (3 2 < fi\ and 



^WO-M > J. > + „; 



2 

Likewise, if v 2 < v] then «2 < «i and 



-u;-u 2 . 



^Kk\ 2 > Jt > (>(a 4 ) + v i 



2 

At any given point one of the above happens. Therefore (14.181) and (14.191) hold with 

fe = ^ (1) («3), 

g e = mm(t/j ( ' (a 4 ) + - , ; (&) + - ) 

and, 

h t = max(^ (2) (a 4 ),^ (2) (/34)). 



Note that for a given e, f e ,g e and fo e are bounded functions. Also for all % = 1,2,3,4, 
-> as e -> and for all i = 1, 2, 3, ft as e ->■ 0. So fl420|) holds. □ 

The Lemma below, is the equivalent of Lemma 14.31 

Lemma 4.8. Let D be a bounded C 2 domain in dimension d > 4. Assume (A2) and that 
di = 0. For x<ED,AcN,3C<oo,e >0, such that Ve < e and y £ 

(4.21) _£|_ £cEA , D(y ,,^ 



Proof of Lemma: As before, one proceeds in two steps. The Step 1 proof, in Lemma 
14.31 follows verbatim. In Step 2, use (14.18j) instead of (12. 20f) to get 

(4.22) v l /(y) < E y {vl'\B Ti )) + E y ( f" h e (B t )v l e (B t )v 2 (B t ) dt). 

Jo 

It is easy to see that the analysis used in obtaining (14.111) . or from [29], will imply that 



for y £ D 8 ' e 



(4.23) AAA ^ Cie A + K *(V> z ^ ■ 
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The induction hypothesis and the fact that h e is bounded will imply 

(4-24) y < c 2 [K x (y, z x + K x (y, z 2 . 

v^yxjv* [x) 

So we have proved the lemma. □ 

We are now ready to state and prove Theorem 14.11 assuming (A2) and not (Al), in the 
case n = 2. 

Theorem 4.3. Let n = 2, and let D be a bounded C 2 domain in dimension d > 4. 
Assume (A2) and that a\ = 0. Then for x,y G D, 

(4.25) lim||| = J^(y^)/ ri = l,2 

(4-26) lim n ^ ( ^ = ^ 2 \0)U {K°(; Zl )K°(;Z 2 )) (y), 

e ~* U LLieA V e\ X ) 

Proof : The proof of (I4.25[) is as in Step 1 of the proof of Theorem 14.11 It essentially 
follows verbatim. Lemma T4.4I was used but that does not require (Al) or (A2). The proof 
of (14.261) follows as in the induction step of the proof of Theorem 14.11 The ingredients for 
identifying the limit and application of dominated convergence are available immediately 
from (jHH]), QUgj) , ( [ODD , <SM, ( M) and <K2T\f respectively. □ 

Remark 4.1. From the above it seems likely that the same idea could work for n > 3. 
To do so would require a suitable appeal to Taylor's Theorem (as in the proof of ft4.18\ )). 
We have not succeeded in carrying this out, and so have presented the proof given earlier, 
under the stronger assumption (Al). 

4.3. Branching backbone for the limiting process, ip analytic. The analysis of 
Section E] can also be carried out for the limiting conditioned process obtained above. We 
will simply state the conclusions one obtains, without repeating the derivation. 

Recall that 

(4.27) M» = Yl \{^\K^), 

where the are given in Theorem 14.11 For e J-^, we let 

1 

be the limiting measure arising in Theorem 14.21 Then the following statements hold, 
under the conditions of that result. 

• is a martingale with respect to J-^,, so the associated Girsanov transform 
defines a consistent probability measure on V^J-fc. In other words, in Dynkin's 
terminology, is an X-harmonic function of X, and is its X-transform. 
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This probability can equivalently be described in terms of a branching backbone 
throwing off mass, that is, as a superprocess with immigration along a random set 
obtained from a branching tree of particles. 

The backbone starts with a single particle, located initially at y. It performs a 
i^-transform of Brownian motion, which dies somewhere in the interior of D. Say 
it dies at y. 

A random partition £ of N is chosen, so given y, the probability that E = a is 



V N (y) 

where a e V(N), \a\ > 2. Here 



v N {y)= KW\^y)\[K*{y)- 

cr&V{N),\a\>2 A€a 

For each A G S, a particle is born at y which proceeds to carry out a i^-transform 
of Brownian motion. If A = {zi}, this particle survives to exit D at z%. If \A\ > 2 
then the particle dies in the interior of D and is replaced by a random number of 
children, labeled by a random partition of A in the manner described above. This 
process repeats until all partitions consist of singletons. This process produces 
a branching tree of particles, with precisely n leaves, corresponding to particles 
exiting D at the n points of N. 

Mass is created/immigrated at points of D where the backbone branches. Given 
that j particles are born because of a branch at y, the mass created is a random 
variable R > whose conditional law ^(dr) is 



r J ir(dr) 
h 1 — S r 2 n{dr) 



J>3 



2a 2 +/ u °r 2 7r(dr)' 3 ~ 

Mass is also created continuously along the backbone according to a Levy process, 
with Levy exponent 

77(A) = tp'(X) - ^'(0) = 2a 2 A + / r(l - e" Ar ) ?r(dr). 

Jo 

Once created, the mass evolves as the (unconditioned) ^-super Brownian motion. 
In other words, with excursion law N.. 



As remarked earlier, note that unlike (13.11) . there is no factorial factor in (14.271) . This is 
simply because of the indexing scheme for the backbones. In Section [3] we ordered the 
n points and then labeled them at random. While now we label using the natural order 
given by the partition. In that sense, the current indexing scheme parallels that of [29J. 

The argument is an induction, based on the Palm formula (Lemma 12.31) . as in Theorem 
13.11 But note that each of the local characteristics of the backbone and mass evolution 
given above are consistent with sending e — > in the characteristics of Section [3] 
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4.4. Explosion of Mass as e — > 0. The explosion effect we wish to understand is already 
present when n = 2, so we will focus mainly on that case. To reiterate the representation 
for analytic tp in this special case, write Kl = K°{-, Zi ). Then the process of interest is 
the Girsanov transform by the martingale 

Ml = (X k ,Kl)(X k ,K*)+m 2 (X k ,U(KlK 2 x )), 

where 



m 2 = — = a 2 + - I r vr(dr). 



oo 

2 



The backbone has a single branch, and follows a [/(i^i^ 2 )-transform till it dies, where- 
upon two particles are born, one doing a ^-transform, and the other doing a K 2 - 
transform. Mass is created at the branch point, according to the law /i 2 , and continuously 
along the backbone, according to the Levy exponent rj. The mass then evolves as an 
unconditioned super-Brownian motion. 

Consider the stable branching function 

POD 

^(A) = / r- (/3+2) (e- Ar - 1 + Ar) dr = C/3 A 1+/3 



o 



(for cp chosen appropriately), which satisfies (12 .4p for < (3 < 1. (Al) fails for ifjp, but it 
does apply to 

POO 

tjP B {\)= / r-^ +2) e^ r (e~ Xr -1 + Ar) dr, 



o 



and clearly ipj — > t^p as 7 \, 0. Thus our construction applies to ipZ, giving an X-transform 
M 2 / 7 with density 

^ = 1 , f 2, 7 = gWVg) +ml{X k ,U(KlKD) 
k Kp\y) k " rnlU{KlKl){y) 

where m] = | / °° r _/3 e~ 7r rfr — > 00 as 7 4- 0. Thus 

which does not represent a valid Girsanov transform, because is not a martingale in 
k under N y . It is not a surprise that something goes wrong here, because expectations of 
terms like (X k , Kl)(X k , K 2 ) should blow up as 7 I - after all, stable random variables 
don't have finite second moments. Our original motivation for carrying out the analysis 
of this paper was understanding precisely what goes wrong in the stable case. In other 
words, of understanding how the singularity arises. 

The problem is not the backbone, since the description of the backbone does not even 
depend on 7. Nor is the problem the continuous mass creation or the subsequent evolution 
of mass, since those approach the corresponding mechanisms for the ^-super Brownian 
motion. The problem is precisely the creation mechanism /i 2,7 (cir) of mass at the branch 
point, since its density is proportional to r~^e _7r , which fails to be tight when 7 4 0. In 
other words, as 7 4 0, the mass born at the branch point blows up. 
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While this heuristic analysis is sufficient to explain the singularity, one can use it to 
rigorously explain the change of measure by the M k . Let ( be the time the original 
particle in the backbone dies. Let Tk be the lesser of the time it dies and the time it exits 
D k . 

Theorem 4.4. M k is a supermartingale. Let <p > 0. Then fory G Dk, 
N y (M k exp-(X k , 0» = limM^exp -(X\ 0), r, < C) 



Proof. Just follow the argument of Theorem 13.11 In the notation of that result, the event 
{ T k < C} is exactly the same as {T k rj ({1,2})}, and the probabilities of such events 
entered into the proof of Theorem 13.11 □ 



The interpretation of the supermartingale property is that we lose absolute continuity 
between the two measures in the limit, with M 2,7 increasingly concentrating mass near a 
set of Ny-measure 0. Namely the set where mass explodes. 

A similar analysis can be carried out in the case n > 2. One shows by induction that as 
7^0, one has b{j,0,y) = (-l)'^(O) ~ 4 j7 1+ ^ , and K A ~ cf^~\ A \U^[ keA Kl) 
for \A\ > 2. It follows that asymptotically there is a single branch point with \A\ branches, 
at which the backbone changes from a single K A transform, to \A\ transforms by K k , 
k G A. Moreover, the mass created at this branch point blows up as 7 -4- 0. For 
example, in the case n = 3, a = {{1}, {2}, {3}} has probability proportional to 7 /3_2 , 
while a = {{1,2}, 3} has probability proportional to 7 /3_1 7 /3_1 = ^ 2 I 3 ~ 2 which is of lower 
order of magnitude. 



4.5. Other orders of limits. Still in the analytic case, with n — 2, the same arguments 
would have handled a slightly more general conditioning. Namely, condition on the exit 
measure charging both and A^, where ei > and 62 > 0. First send ei — > and 
then €2 — > 0. It is not hard to check that this gives precisely the same limiting object 
as before. So if we apply this procedure to the ^-super Brownian motion, and then 
let 7 J, 0, the mass should still blow up, but the backbone should take the form of a 
[/(iT^.K^)-transforrn, splitting into a ^-transform and a ^-transform. 

With this modification, one should be able to treat a different approach to the stable 
case than given above. Namely, start out with the actual ^-super Brownian motion. 
Condition its exit measure to charge A*J and A^. Then send e± — > followed by e 2 — > 0. 
The following informal calculation suggests how the backbone should behave in the limit. 
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For % — 1, 2 set 

u* = v* = < = < = M.(X D (Af;) > 0) 
n 12 = < 2 £2 = N.(X D (A £ 2 i 1 U A e z 2 2 ) > 0) 

w 12 = v \l €2 = u l + u 2 - u 12 = K(X D (A z e l) > 0,X D (A Z2 2 ) > 0). 

Then 

^At; 12 = l -A[u l +u 2 - u 12 } = M^ 1 ) + Mu 2 ) ' Mu 12 ) 

= Mv 1 ) + IMv 2 ) - M"" 1 + v2 )} + [Mv 1 + - Mv 1 + v 2 - v 12 )} 
~ ^(V] 1+/3 - [i + 0[« a ]V + [i + fll* 1 + v'lV 2 )- 

Fix x and rescale, letting v 1 = v 1 (-)/v 1 (x), v 2 = v 2 (-)/v 2 (x), v 12 = v 12 (-)/v 1 (x)v 2 (x) 13 . If 
all these quantities remain bounded, then sending e\ — > should give 

l -/\v 12 ~ (1 + /3) C/3 (V]V 2 - [t; 2 ]^ 1 ). 

So sending ei — >■ and then 62 — ► should give 

i A D 12 ~ -(1 + Pfo^v 1 . 

Though we again expect this conditioning to have a mass which blows up in the limit, the 
above suggests that the backbones should still converge weakly. But this time the limit 
should be a C/([i^ 2 ]^^)-transform, splitting into a ^-transform and a .^-transform. 

Call the first part of the backbone the trunk. Taking limits in the other order should 
produce a trunk that is a ^([i^]^i^ 2 )-transform. So taking limits in which both e± — > 
and €2 — > in a coordinated way should give trunks that are transforms by 

u([KlK 2 n9Kl + (l-9)K 2 ]^), 

for arbitrary < 9 < 1. 

We will not try to make the above informal argument rigorous. But we originally found 
it puzzling. Not so much because uniqueness of limits breaks down in the case of the 
^a-super Brownian motion. But rather because the U(K],K 2 ) trunk we obtained earlier 
does not appear among the possible limits when taken this other way. This suggests that 
there should be a way of capturing a broader class of limits, encompassing both types 
obtained above. Or at least, of obtaining both types of limits by a common procedure. 
Carrying this out rigorously, in the current context, seems more work than it is perhaps 
worth, and we do not claim to have done so in full detail. But in the next section, we will 
find a simpler setting, in which this can be done more more easily. 



BLOWUP AND CONDITIONINGS OF V-SUPER BROWNIAN EXIT MEASURES 



41 



4.6. A class of martingales, with n — 2. In this section, we focus on a technically 
simpler collection of martingales, exhibiting some of the same behaviour as found above. 
We only present the heuristic idea and do not present detailed proofs in this section. The 
heuristic development can be made rigorous but we have chosen not to include this, as our 
aim in this section is to illustrate how a spectrum of limits can be obtained via various 
limiting procedures.. 

Let if; satisfy (Al). Let it, f, g, v > be bounded solutions on D to the following equations: 



(4.28) 





1 A 

-Ait 

2 


= ip(u) 


2 J 




= 


% 

2 y 


- *l>\u)g 


= 


1 . 

-Av 

2 


— ip'{u)v 


= -r 



as well as v = on dD (so v is an L^/ ou -potential). Call M.2 the set of (u,f,g,v) so 
obtained. 



Under ML,,, for every (it, /, g, v) £ M. 2 define: 



M k (u, f, g, v) = e-<^> \(X k , v) + (X k , f) (X k , g) 



It can be verified that Mfe(it, /, g,v) is a T k martingale. Moreover, we can realize the 
Girsanov transform 

1 

v(y) 



can be described as 



by such martingales, as follows: Let (it, f,g,v) £ A^2- Then 
follows: Start a i>-transform of the L^i ou process till it dies in D. At this point, start two 
paths that run to 3D, respectively an /-transform and a (yf-transform of the L^i ou process. 
Mass is created at the branch point, according to the law /x 2 , and continuously along the 
backbone, according to the Levy exponent 77. The mass then evolves as an unconditioned 
super-Brownian motion. 



In this simplified context, the analogue to the question considered in previous sections 
is the following: Let (u n , f n , g n ,v n ) £ M.2, with each term — > 0. Find the limit points, 
either of M^ nJn ' 9n ' Vn) , or of the backbone. 

When ip satisfies (Al), there is a systematic answer to this question. For a fixed x, we 
renormalize by constants a n and b n , so that 
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converge to non-zero values. Then f n , g n , and v n = -^r- all have non-zero limits f,g,v 
which satisfy 

±a/ = ±a</ = o, v = r(o)u(fg)- 

Moreover M^ n ' Sn ' 9n ' Vn) converges in the weak topology to M { y' f,g,v)) . 

In the case of ipp, the actual Girsanov transforms yi y Un >f n > 9n > Vn "> degenerate as before, 
because mass blows up. We look at the limits of the backbones instead. 

Let a n , b n , c n all be sequences of positive reals, that all — > 0. Assume that a n b n = o(c n ~P). 
Let <pf,<j)g,<f) u be smooth functions on dD, that are bounded away from 0. Define 
(u n , f n , g n ,v n ) to be the solutions to (I4.28p . using boundary conditions a n (f>f for f n , b n 4> g 
for g n , and c n <f)g for u n . Set d n = + /3)a n b n c n ^ and 

_ 7 fn — 9n — _ 

u n i in j Qn T i 

C n d n O n Cl n 

(14281) gives that 

^Av n - (1 + p)u£v n = -0(1 + f3)u-^a n bj n g n , 

so 

^Av n - (1 + P)vfc) n = -U- {1 -^fn9n. 

Now taking limits we have u n — >• it, / n — > /, g n — > g, and v n — > v where u, f, and 
g are the harmonic functions with boundary values (p u , (fif, (f) g , and v = U(u~^ l ~^ f g). 
The backbones converge weakly to a ^-transform branching into an /-transform and a 
^-transform. 

We conclude this section with two of examples analogous to the limits obtained in the 
previous sections. 

Example 4.1. Take <p u = 1. The trunk is a U (f g) -transform, as in Section \4^\ 
Example 4.2. Take <p u = <pf. The trunk is a U (f ^ g) -transform, as in Section \47h*\ 

In other words, this gives a general class of limiting objects, which encompasses examples 
analogous to both types of limits obtained earlier. In that sense, the model of this section 
interpolates between these examples, and helps explain the variety of limits obtained. 

References 

[1] Abraham, R. and Le Gall, J.F. (1994) Sur la mesure de sortie du super mouvement brownien. 

Probabability Theory and Related Fields 99, 251-275. 
[2] Azema, J. (1979) Theorie generate des processus et retournement du temps. Ann. Sci. Ecole Norm. 

Sup. 6, 459-519. 



BLOWUP AND CONDITIONINGS OF V-SUPER BROWNIAN EXIT MEASURES 



43 



[3] Banuelos, R., Bass, R., and Burdzy, K. (1991) Holder domains and the boundary Harnack principle. 

Duke Math. J. 64, no. 1, 195-200. 
[4] Chung, K. L. and Zhao, Z. (1995) From Brownian motion to Schrddingcr's equation. Springer- Verlag, 

Berlin. 

[5] Cranston, M., Fabcs, E., and Zhao, Z. (1988) Conditional gauge and potential theory for the 
Schrddingcr operator. Trans. Amer. Math. Soc. 3307, 171-194. 

[6] Dawson, D. A. and Perkins, E. A. (1991) Historical Processes. Memoirs of the American Mathemat- 
ical Society, No. 454. 

[7] Dawson, D. A. (1991) Measure-valued Markov Processes. Ecole d'Ete de Probabilites de Saint-Flour 

XXI, Lecture Notes in Mathematics 1541, Springer- Verlag New York. 
[8] Dynkin, E.B. (1991) An Introduction to Branching Measure-Valued Processes. CRM Monographs 

series 6: Amer. Math. Soc, Providence, R. I. 
[9] Dynkin, E. B. (1991) A Probabilistic Approach to One Class of Nonlinear Differential Equations. 

Probability Theory and Related Fields 90, 1-36 
[10] Dynkin, E. B. (2002) Diffusions, Superdiffusions and Partial Differential Equations, Amer. Math. 

Soc, Providence, R.I. 

[11] Dynkin, E. B. (2004) Harmonic functions and exit boundary of supcrdiffusion, Journal of functional 
analysis 206, 33-68. 

[12] Dynkin, E. B. (2004) Superdiffusions and positive solutions of nonlinear partial differential equations. 
University Lecture Series 34, Amer. Math. Soc, Providence, R. I. 

[13] Englander, J. and Pinsky, R. (1999) On the construction and support properties of measure-valued 
diffusions ouflCR with spatially dependent branching. Ann. Probab. 27, no. 2, 684-730. 

[14] Ethcridgc, A.M. (1993) Conditioned superprocesses and a scmilincar heat equation, in "Seminar on 
Stochastic Processes 1992" , Birkhauser, Boston, 9199. 

[15] Etheridge, A. M. and Williams, D. R. E. (2003) A decomposition of the (1 + f3)-superprocess condi- 
tioned on survival. Proc. Roy. Soc. Edinburgh Sect. A 133, 829-847. 

[16] Evans, S.N. (1993) Two representations of a conditioned supcrprocess. Proc. Roy. Soc. Edinburgh 
Sect. A 123, 959-971. 

[17] Evans, S.N. and Perkins, E.A. (1990) Measure-valued Markov branching processes conditioned on 

non-extinction. Israel J. Math, 71: 329-337. 
[18] Grey, D. R. (1974) Asymptotic behaviour of continuous time, continuous state-space branching 

processes. J. Appl. Probab. 11, 669-677 
[19] Kyprianou, A.E., Liu, R.-L., Murillo-Salas, A., and Ren, Y.-X. (2010) Supercritical super-Brownian 

motion with a general branching mechanism and travelling waves. Preprint. 
[20] Le Gall, J.F. (1994) Hitting probabilities and potential theory for the Brownian path-valued process. 

Ann. Inst. Fourier (Grenoble) 44, 277-7306. 
[21] Le Gall, J.F. (1996) Superprocesses, Brownian snakes and partial differential equations. Prcpubli- 

cation, Lab. de Probabilites, Univ. Paris VI 337. 
[22] Meir, A. and Moon, J.W. (1978) On the altitude of nodes in random trees. Canad. J. Math. 30, no. 

5, 997-1015. 

[23] Moras, M. (2009) Conditioned super Brownian motion in Denjoy domains and strips. Ph. D. Thesis, 
York University. 

[24] Overbeck, L. (1993) Conditioned super-Brownian motion. Probab. Theory Relat. Fields, 96, 545-570. 

[25] Overbeck, L. (1994) Pathwise construction of additive H-transforms of super-Brownian motion. 
Probab. Theory Relat. Fields 100, 429437. 

[26] Perkins, E. A. (2002) Dawson-Watanabe superprocesses and measure-valued diffusions. Lecture 
Notes in Mathematics, 1781. Springer- Verlag, Berlin, (125-329). 

[27] Roclly-Coppolctta, S. and Rouault, A. (1989) Processus de Dawson-Watanabe conditionne par le 
futur lointain. C.R. Acad. Sci. Paris, Scr. 1 309, 867-872. 

[28] Salisbury, T.S. and Sezer, A.D. (2010) Conditioning super-Brownian motion on its boundary statis- 
tics; fragmentation and a class of weakly extreme X-harmonic functions. Preprint. 



44 



SIVA R. ATHREYA AND THOMAS S. SALISBURY 



[29] Salisbury, T. S. and Verzani, J. (1999) On the conditioned exit measures of super Brownian motion. 

Probability Theory and Related Fields 115, 237-285. 
[30] Salisbury, T. S. and Verzani, J. (2000) Non-degenerate conditionings of the exit measures of supcr- 

Brownian motion. Stochastic Processes and Their Applications 87, 25-52. 
[31] Serlct, L. (1996) The occupation measure of supcr-Brownian motion conditioned on noncxtinction. 

J. Thcorct. Probab. 9, 561-578. 
[32] Veron, L. (1996) Singularities of solutions of second order quasilinear equations. Pitman Research 

Notes in Mathematical Series 353, Addison Wesley Longman Ltd, U.K. 
[33] Verzani, J. (2008) (A class of extreme X-harmonic functions.) Commun. Stoch. Anal. 2, 335-355. 



Siva R. Athreya, 8th, Mile Mysore Road,, Indian Statistical Institute, Bangalore 560059 
India. 

E-mail address: athreyaOisibang. ac . in 

Thomas S. Salisbury, Dept. of Math & Stat., York University, Toronto, Ontario, Canada 
M3J 1P3 

E-mail address: salt@yorku.ca 



